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We have studied the CP violation discovery reach in neutrino oscillation experiment with su¬ 
perbeam, neutrino factory and monoenergetic neutrino beam from electron capture process. For 
NSI satisfying model-dependent bound for shorter baselines (like CERN-Frejus set-up ) there is 
insignificant effect of NSI on the the discovery reach of CP violation due to 5. Particularly, for 
superbeam and neutrino factory we have also considered relatively longer baselines for which there 
could be significant NSI effects on CP violation discovery reach for higher allowed values of NSI. 

For monoenergetic beam only shorter baselines are considered to study CP violation with different 
nuclei as neutrino sources. Interestingly for non-standard interactions - Seii and Ser of neutrinos with 
matter during propagation in longer baselines in superbeam, there is possibility of better discovery 
reach of CP violation than that with only Standard Model interactions of neutrinos with matter. 

For complex NSI we have shown the CP violation discovery reach in the plane of Dirac phase S 
and NSI phase (pij. The CP violation due to some values of S remain unobservable with present 
and near future experimental facilities in superbeam and neutrino factory . However, in presence 
of some ranges of off-diagonal NSI phase values there are some possibilities of discovering total CP 
violation for any 6cp value even at 5cr conhdence level for neutrino factory. Our analysis indicates 
that for some values of NSI phases total CP violation may not be at all observable for any values 
of S. Combination of shorter and longer baselines could indicate in some cases the presence of NSI. 
However, in general for NSIs < 1 the CP violation discovery reach is better in neutrino factory 
set-ups. Using a neutrino beam from electron capture process for nuclei 50 °Sn and ^®^Yb, we have 
shown the discovery reach of CP violation in neutrino oscillation experiment. Particularly for so^Sn 
nuclei CP violation could be found for about 51% of the possible 5 values for a baseline of 130 km 
with boost factor 7 = 500. The nuclei ^®^Yb is although technically more feasible for the production 
of mono-energetic beam, but is found to be not suitable in obtaining good discovery reach of CP 
violation. 
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I. INTRODUCTION 


Long back CP violation has been found in the quark sector of the Standard Model of Particle Physics. But 
so far there is no evidence of CP violation in the leptonic sector. One way to search for such CP violation 
is through neutrino oscillation experiments in which one flavor of neutrino could oscillate to other flavors of 
neutrinos. Neutrino oscillation probability depends on various oscillation parameters present in the 3x3 
mixing matrix - the PMNS matrix [1-3] and the neutrino mass squared differences. Two of the three angles - 
012 and 023 present in PMNS matrix have been known with certain accuracy for some time. Recently several 
experiments like Double CHOOZ, Daya Bay, RENO and T2K collaboration [4-7] found non-zero value of 
sin^ 2013 corresponding to the third mixing angle of PMNS matrix with a global level of significance which 
is well above conventional 5 ct discovery threshold. As three angles are non-zero there could be non-zero CP 
violating phase 5 in the PMNS matrix and as such there could be CP violation in the leptonic sector. The 
mass squared differences - [Aml^] and is known (where Amfj = mf—mj) but the sign of Aml^ and as 

such the hierarchy (whether it is normal (NH) or inverted (IH)) of neutrino masses is still unknown. Various 
neutrino oscillation experiments like superbeam, neutrino factory, beta beams and reactor experiments are 
focussing on determining these unknown parameters corresponding to neutrino oscillations [8-19]. 

There could be various kind of non-standard interactions of neutrinos with matter. This could have some 
effect on finding CP violation in neutrino oscillation experiments. Such interactions could be the non¬ 
standard interactions of neutrinos with matter fermions {u, d and e) during propagation of neutrinos only. 
This could affect oscillations of different flavors of neutrinos as sub-leading effect. We have discussed it in 
further detail in the next section. There could be other different kind of interactions beyond Standard model 
leading to non-unitarity of 3 x 3 PMNS neutrino mixing matrix [20]. Considering non-standard interactions 
of neutrinos at the source and the detector in neutrino oscillation experiments also lead to such possibility. 
However, such NSI at the source and detector have highly stringent constraints [21] and as such the effect 
on neutrino oscillation is expected to be lesser affected than that due to NSI in matter during propagation 
of neutrinos. We shall not consider NSI at the source and detector in this work. 

To study CP violation in neutrino oscillation experiments we have considered superbeam, neutrino factory 
and monoenergetic neutrino beam from electron capture process. There are some earlier studies for short and 
long baselines for standard interactions [18, 19, 22-26] for neutrinos coming from superbeam. We consider 
neutrino superbeam (which mainly contains and ) coming from CERN travelling a baseline length of 
2300 km to Pyhasalmi (Finland) and another baseline of 130 km to Frejus (France). We have also considered 
superbeam coming from Tokai to Hyper-Kamiokande travelling a baseline of 295 km. We do a comparative 
study in the discovery potentials of the CP violating phase 6 for these three baselines in presence of both 
standard and non-standard interactions. Main emphasis of our work is on exploring CP violation due to 5 
in the light of recent experimental findings of large value of 0i3 and how NSI could affect the discovery reach 
of such CP violation. Also there are some studies on CP violation in presence of non-standard interactions 
[27, 28] before Daya Bay experiment when 0i3 was not so precisely known. 

The discovery of CP violation for different parent muon energy in neutrino factory and different length of 
baselines have been studied earlier [29, 30]. In the analysis done in [29] three different kinds of magnetized 
detectors: a Totally-Active Scintillator(TASD) and two kinds of Liquid Argon detectors have been considered 
for baselines ranging from 1000 km to 4000 km and parent muon energies from 4 GeV to 25 GeV, whereas in 
the analysis done in [30] MIND (Toroidal magnetized iron neutrino detector) detector has been considered. 
As recent reactor neutrino experiments indicates large value of 0i3 it is important to study the discovery 
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potential of CP violation in neutrino oscillation experiments in low energy neutrino factory [31]. There 
are some studies on the performance of low energy neutrino factories [15] in the context of standard model 
interactions of neutrinos with matter [29, 31-38] and also in presence of NSI [39-43]. In the works related 
to NSI [39-43] longer baselines like 4000 and 7500 kms has been considered to have larger NSI effect on 
neutrino oscillation experiments so as to find better NSI sensitivity and in that context CP violation has 
also been studied. 


Our work is complementary to these earlier analysis in the sense that we focus here on the better discovery 
reach of CP violation and as such we consider relatively shorter baselines where the matter effect and as 
such NSI effect will be lesser in the neutrino oscillation experiments. However, in the context of such 
experimental set-up in neutrino factory which is better optimized for CP violation discovery (considering 
only SM interactions of neutrinos with matter) we have also addressed the question oi CP violation coming 
due to NSI phases. We have considered lower parent muon energy within 10 GeV which is found to be more 
favourable for the discovery of CP violation. The optimization for CP violation discovery in the context of 
SM interactions has been done earlier [29, 30]. However, in this work we re-analysed this optimization with 
the updated detector characteristics for the MIND detector [44]. Based on this optimization we have chosen 
a few relatively shorter baselines (730, 1290 and 1500 kms) and low parent muon energy of 10 GeV. 


In recent years oscillation experiment using a neutrino beam with neutrinos emitted from an electron 
capture process is proposed [16, 45-51]. Such beam can be produced using an accelerated nuclei that decay 
by electron capture. In this process an electron is captured by a proton releasing a neutron and an electron 
neutrino. So the beam is purely of one flavor. In the rest frame of the mother nuclei the electron neutrino 
that is released from such process, has a definite energy Q. Since the idea of using a neutrino beam emitted 
from an electron capture process is based on the acceleration and storage of radioactive isotopes that decays 
to daughter nuclei, one may get the suitable neutrino energy by accelerating the mother nuclei with suitable 
Lorentz boost factor 7 . One can control the neutrino energy by choosing the appropriate Lorentz boost 
factor as the energy that has been boosted by an appropriate boost factor towards the detector is given as 
E = 2'yQ. Hence for certain mother nuclei to get the required neutrino energy the boost factor have to 
be chosen appropriately with respect to Q. Due to the almost monoenergetic nature of such beam one can 
appropriately choose the neutrino energy for which the probability of oscillation could be large and sensitive 
to certain unknown neutrino oscillation parameters. 


In this work we consider such a flavor pure electron neutrino beam emitted from electron capture process 
for a suitable 7 value where the beam is targeted towards a Water Cherenkov detector and perform numerical 
simulation to study the discovery reach of CP violation in oscillation experiments. 


In section H we have mentioned various model dependent and independent bounds of NSIs and have 
discussed the neutrino oscillation probability in presence of NSI. In section HI, we have discussed the pro¬ 
cedure of numerical simulation and in its different subsections related to superbeam, neutrino factory and 
monoenergetic neutrino beam we have discussed various experimental set-ups and results of our analysis on 
discovery of CP violation. In section IV we have concluded with some remarks on the possibility of knowing 
the presence of NSIs in the context of superbeam and neutrino factory. For monoenergetic beam we have 
discussed the suitable nuclei for discovering CP violation and related technical constraints. 
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II. NSI AND ITS EFFECT ON NEUTRINO OSCILLATION PROBABILITIES 


We consider the non-standard interactions of neutrinos which could be outcome of effective theory at low 
energy after integrating out the heavy mediator fields at the energy scale of neutrino oscillation experiments. 
Apart from Standard Model (SM) Lagrangian density we consider the following non-standard fermion- 
neutrino interaction in matter defined by the Lagrangian: 


£^si = -2V2GFefj;[h^Pf][DprL,.^] ( 1 ) 

where P G {L,R), L = \ R = ^ , / = e,u,d and are termed as non-standard interactions 

(NSIs) parameters signifying the deviation from SM interactions. These are non-renormalizable as well as 
not gauge invariant and are dimension-6 operators after heavy fields are integrated out [21]. Although at 
low energy NSI look like this but at high energy scale where actually such interactions originate there they 
have different form. These NSI parameters can be reduced to the effective parameters and can be written 
as: 


— 


E 


a/3 „ 


( 2 ) 


where n/ and Ue are the fermion and the electron number density respectively in matter. As these NSIs 
modify the interactions with matter from the Standard Model interactions the effective mass matrix for the 
neutrinos are changed and as such there will be change in the oscillation probability of different flavor of 
neutrinos. Although NSIs could be present at the source of neutrinos, during the propagation of neutrinos 
and also during detection of neutrinos [52] but as those effects are expected to be smaller at the source and 
detector due to their stringent constraints [21, 53], we consider the NSI effect during the propagation of 
neutrinos only. 

Model dependent [21, 54-66] and independent [53, 67, 68] bounds are obtained for these matter NSI 
parameters and are shown in the following table. In obtaining model dependent bounds on matter NSI the 
experiments with neutrinos and charged leptons - LSND, CHARM, CHARM-II, NuTeV and also LEP-II have 
been considered. Bounds coming from loop effect have been used for model dependent bounds. However, 
model independent bounds are less stringent and could be larger than the model dependent bounds by several 
orders and have been discussed in [21, 53]. Considering recent results from experiments in IceCube-79 and 


NSI 

Model dependent 
bound on NSI [Reference [21]] 

Model independent 
Bound on NSI [Reference [53]] 

^ee 

> -0.9; < 0.75 

< 4.2 


^ 3.8 X 10-" 

< 0.33 

\^eT 1 

^ 0.25 

< 3.0 

Sn/i 

> -0.05; < 0.08 

< 0.068 

k/iT 1 

^ 0.25 

< 0.33 

S-TT 

^ 0.4 

< 21 


TABLE I: Strength of Non standard interaction terms used for our Analysis. 
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DeepCore more stringent bound on \e^r \ and Ett have been obtained in [ 68 ]. In section IV, we shall 
consider both model dependent and independent allowed range of values of different NSIs as shown in the 
table above for earth like matter while showing discovery reach for CP violation in presence of NSI. 


A. Ve —> Vfi oscillation probabilities with NSI 


The flavor eigenstates Va is related to mass eigenstates of nentrinos Vi as 

^OL — ^ ^ U ai^i 5 ^ — 1) ^5 Hi (3) 

i 

in vacuum where 17 is 3 x 3 unitary matrix parametrized by three mixing angles di 2 , O 23 and ^13 and by one 
CP violating phase 5 for Dirac neutrino as given in [2, 3]. However, for Majorana neutrino apart from 5 
there are two more CP violating phases which do not play role in neutrino oscillation [ 2 ]. 

In the flavor basis the total Hamiltonian consisting both standard (Hsm) and non-standard interactions 
(Hnsi) of neutrinos interacting with matter during propagation can be written as: 


H = Hsm + H^si 


where 


Hsm 


Awji 

2E 


/o 0 o\ 

H 0 a 0 1/1^ 

. VO 0 ij 


/a 0 oV 
000 , 

\o 0 oy 


77nsi 


2E 


A 


(e 

c-ee 



In equations (5) and ( 6 ) 


2E'/2GFne Am^i 


Am'fj = m} — rn?j 


( 4 ) 


( 5 ) 

( 6 ) 


( 7 ) 


where rrii is the mass of the z-th neutrino, A corresponds to the interaction of neutrinos with matter in SM 
and Gf is the Fermi constant. See, ^^0 £tt correspond to the non-standard interactions 

(NSIs) of neutrinos with matter. In equation ( 6 ), ( * ) denotes complex conjugation. The NSIs - Ee^, £er 
and e^t could be complex. Later on, in the expressions of probability of oscillation we have expressed these 
NSIs as Eij = \Eij\P‘^'C In our numerical analysis we have considered the NSIs - Ee^, Eer and e^r as both 
real as well as complex. 
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For longer baselines following the perturbation method adopted in references [ 69 , 70 ] for the oscillation 
probability upto order o? (considering sin0i3 ~ \foi as follows from recent reactor experiments) and 

the matter effect parameter A in the leading order of perturbation (which happens for longer baselines) and 
NSI parameters Sap of the order of a one obtains [ 71 ] 


Pu 


= P: 


•SM 


P. 


MSI 


( 8 ) 


where 


P^^,^ = 4 sin 


(A l)Aml,L slssls -a+ A)hl3) sin l)AmiiL 


,2 r „2 2 

'°13'’2; 

AE {A-l)‘^ 


AE 
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A{A - cos -^ sm^ sm^ 

Q^gi2gi3'S23 ( {A-l)Aml^L {A-l)Aml^L _ „ . . 2 {A-l)Aml^L 
( 21 - 1)3 \ E 2 E AE 

asi 3 S 2 xi 2 S 2 x 23 /c Am§iL\ , ( 21 - 1 )Ato|iL . AAm'^iL 

A{A- 1 ) \ \ AE J AE AE 


( 9 ) 


P, 


NSI 


4|a2|s2x23Si3 . AAm'^iL . (2l-l)Am§iL 


2l(2l - 1 ) ■ AE 
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(A - 1)2 
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(|a3| + 2 cos ((5 + ^03)813) 


( 5 - 


Arri^iL 
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where 


ai = As 

|a 2 |e *'^“2 = A[ e*'^'^|ee^|c23 - |eer|s23 


lasle* “3 ^6x1^23 + ^^\^efi\s23 

|a4|e*'^“'‘ = — 2|e^T-|s23 + {Sfj,fj, — eTr)c23S23^ 

rt5 — ^f^^rTC23 “t“ l^/xr | COS (/)^'7-52x 23 


( 11 ) 


and 


(j}a2 = tan“^ 
(/)a3 = tan“^ 
(j)ai = tan“^ 


|£eM|c 23 Sm(/)e^ - [ger | ^23 sin (?^eT 

I C23 COS ^e/i] l^-er I COS S23 _ 

l^e/j-1 ^23 sin. (pefi + |£er|c23 sin (/) er 1 

I'^er |C23 COS(/) er “1“ l^e/j, I cos (j)etiS23\ 

|£;/t| sin(^^ 7 -) 

|£^tT|c2x23 C0s((5^>y3x) “1 “ £tt)c23'^23 


( 12 ) 


where Sij = sinOij, Cij = cos Oij, S2xij = sin2dij, C2xij = cos20y . The oscillation channel —>■ is 
important for neutrino factory. For superbeam, —t Ve oscillation channel is important which can be 
obtained by the following transformation 

Pa/3iA, S, <j)ij) = P/3a{A,-6, -(pij). (13) 

The oscillation probabilities for antineutrinos can be obtained from the oscillation probabilities given for 
neutrinos above by using the following relation: 

P^-^{S, A, = P^^{-S, -A, (14) 


To estimate the order of magnitude of 6 dependent and 6 independent but matter dependent ( i.e., A 
dependent) part in the above two oscillation probability, following results of reactor experiments we shall 
consider sin 0 i 3 ~ ^/a. For only SM interactions, (i.e £q ,/3 —>■ 0) from last term of Eq.(9) for oscillation 
probabilities one hnds that the 5 dependence occurs at order of asi 3 ~ for both neutrino oscillation 
and anti-neutrino oscillation probabilities. 

In presence of NSI of order a, in the first two leading terms of Eq.(lO), the <5 dependence occurs at order 
of 02513 ~ 03513 ~ through terms containing 02 and 03 (which are Se^i and Ser dependent). Due to this 
(if we assume all NSI of same order ~ a), these two NSI {Seji. and £er) in contrast to other NSI could have 
greater effect on oscillation probability in —)■ v^^ or the reverse oscillation channel . 



For shorter baselines the above oscillation probability expression becomes much simpler. Following the 
perturbation method [52] and by considering the standard model matter effect A ^ a for neutrino energy 
E around 1 GeV and sin0i3 ~ ^/a as obtained from reactor oscillation data, the probability of oscillation 
upto order is given by 


where 


P. 


= P. 


, (a) + Pv^^u^ + Pv e—) 


(15) 


Pv^^u. (a) = sin^ 2013 sin^ 023 sin" 


P. 


P.. 


( AmoiL 

sin 013 sin 2023 sin20 
/ Ato|iL^ ^ 


. / AtooiL 


, - AmhL 


^ (a^) = cos^ 023 sin"^ 20i2 


V 4f; 


— 2q; sin 013 sin 0 i 2 sin 623 


Arri^iL . ArnX^L 


2E 


sm ■ 


31J 


2E 


+ 8A sin^ 013 sin^ 623 ( sin' 


Am§iL Am^iL . ArrX-iL 


31^ 


AE 


8 E 


sm ■ 


‘31-* 


2E 


(16) 


One may note that this expression of oscillation probability is little bit different from that presented by 
earlier authors [72] because they have considered the perturbative approach for relatively longer baseline 
and small sin0i3 ~ a. For shorter baseline as presented here, the S dependence appears to be more and 
instead of being at the order of oX it appears at the order as a sin 0 i 3 ~ although the expression 
is same. However at order cX some extra terms appear in comparison to that presented by earlier authors 
although their contribution is relatively smaller being at the order of oX. The matter effect is occurring at 
order oX through term containing parameter A. 


III. NUMERICAL SIMULATION 


There is detailed global analysis of three flavor neutrino oscillation data [73] presented at the Neutrino 
2012 conference. In their analysis the correlation between various oscillation parameters has been taken into 
account. Very recently there is another such global analysis [74] presented at the Neutrino 2014 conference 
which has been considered in our analysis and best fit values of neutrino mixing parameters and their 
respective errors at 3a confidence level are shown in table II. For earth matter density the PREM profile 
[75] has been considered. Also we have considered an error of 2 % on matter density profile. 


Oscillation Parameters 

Central Values 

Error {3a) 

^”^21 py 2 

in-5 ev 

7.45 

7.02 - 8.09 

10-3 ev 

2.457 

2.325 - 2.599 (NH) 
(-2.59)- (-2.307) (IH) 

sin"^ 012 

0.304 

0.270 - 0.344 

sin"^ 023 

0.452 

0.385 - 0.644 

sin"^ 013 

0.0218 

0.0188 - 0.0251 


TABLE II: Central values of the oscillation parameters with errors. 
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The numerical simulation has been done by using GLoBES [76, 77]. GLoBES uses poissonian for 
the oscillation parameters A and for nuisance parameters To implement the systematic errors the ‘pull 
method’ [78] has been used by GLoBES and Xpuii is given as [77] 

xluiiW :=min("x^(A,Ci,--- + ( 1 ’^) 

G.l V CTj, / 

where X^(A,^i, • • • is the usual Poissonian x^ depending on the neutrino oscillation parameters A and 
the nuisance parameter and it has been summed over different energy binning on the theoretical and 
observed event rates. The oscillation parameters A is 

= {(^12, (^13, (^23, Amh, p,ScP,£al3,(l)al3) ; (18) 

The nuisance parameters are the signal and background normalization and the calibration errors. In order 
to implement the error , the x^ (in which all event rates of all bins are included) is minimized over the 
different nuisance parameters independently. 

We have also marginalized over A' where 


A' = (6»i2,6»i 3,6»23, Am2i, Amgi.p) . 


(19) 


The final projected Xf i® given as 


x| = inin (xp„o(A) + priorsiX')) 


( 20 ) 


where 


priors{X') = 


P- P 


/ sin^ 9ij - sin^ 9^^ ^ ^ 




O'sin^ 0i 


Am|i - {Aml^f 


+ 


|Am§i| - |(Am §i 
o-|Am§j 


2 N0| 


( 21 ) 


and sin^ 6 *]' 2 , sin^^Jg, sin^ 6 * 23 , (Am 2 i)°and|(ATO§ 2 )°|) are the central values of the corresponding pa¬ 
rameters. We have marginalized over two different hierarchies of neutrino masses also. 


For taking into account NSI in GLoBES we have followed the method described in [77] and modified 
the source file ”probability.c” in GLoBES appropriately by inserting the NSI’s in the subroutine where 
the hamiltonian for matter interaction is defined. Then we have included the new probability program as 
instructed by the manual of GLoBES. We have not used the approximate analytical expression of oscillation 
probability for NSI of order a in our numerical analysis. 


At present we do not know the hierarchy of neutrino masses and the value of 6. In the test values we have 
considered both possibilities of the hierarchy. However, we shall show the results in Figures only for true 
normal hierarchy. Before we get the knowledge of specific value of 6 it could be possible to know whether 
nature admits CP violation or not which could be possible for a wide range of values of 6. So what we are 
going to analyse is not the measurement of S values but finding out all possible d (true) values which could be 
distinguishable from the CP conserving S values at certain confidence level through oscillation experiments. 
For that it is useful to calculate the CP fraction due to 6 denoted as F{6) which may be defined as the 
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fraction of the total allowed range (0 to 2 tt) for the 5(true) over which CP violation can be discovered 
[29, 40]. While calculating the fraction due to d for the discovery of CP violation, we consider the exclusion 
of all parameter sets with 5 G { 0 , 7 r} which means (5(test) is fixed in its CP conserving values { 0 , 7 r}. For 
no NSIs, the degrees of freedom in analysis is one due to due to <5. We have calculated the minimum of 
those for each (5(true) and have calculated the fraction of the i5(true) over its entire range which satisfies 
> 9(25) at 3(5 )ct confidence level. 


In studying the NSI effect (one at a time) on CP violation we have fixed absolute values of NSI(test) to 
their true values which implies that the absolute values of NSI are assumed to be known from some other 
experiments. In Figures we have shown the effect of real NSI’s on the CP fraction discovery. In this case, 
the degrees of freedom in x^ analysis is one due to 5. The off-diagonal NSI in NSI matrix could be complex 
and the corresponding complex phases could be new sources of CP violation. If the CP violation is only 
due to NSI phase, then the CP fraction due to NSI phases (j)ij (denoted as F{(j)ij)) is calculated in the same 
way as described above for S by replacing S by the corresponding (j>ij. Here the degrees of freedom is one 
due to (j)ij. In Figures we have shown the effect of absolute values of non-diagonal NSI’s (one at a time) on 
F{(j)ij). In order to calculate the the total CP violation discovery due to two phases (pij and 6 we have set 
(j)ij{test) and (5(test) to their respective CP conserving values at {S,(j)ij) € {(0,0), (0, tt), (tt, 0), (tt,tt)} and 
also both of the mass hierarchy in the test values have been taken into account. Here, the degrees of freedom 
is two due to 6 and one of the NSI phases (pij. Then we have calculated the above x^ minimum for each 
true values of (fiij and 6 and have shown the allowed region for discovery in (j)ij-6 space for > 11.83 at 
3 ct confidence level. In all Figures true hierarchy has been considered to be normal. 


So by numerical simulation we are studying the different possibilities of finding CP violation in different 
experimental set-ups due to 6 or NSI phase (fiij or total CP violation due to S and NSI phase (j)ij for certain 
absolute value of NSI (e^ (true)) assumed to be present in nature. If we assume that any allowed values of 
d has equal probability to be true 5 value of nature then lower (higher) CP fraction which we present in the 
result sections will indicate lower (higher) probability of finding CP violation in particular experimental- 
set-up discussed below. One may note that NSI effect and the 5 effect in the hamiltonian H has different 
neutrino energy dependence (see equation Eq.(5) for <5 present in U and Eq.( 6 ) for NSI). Also in the shorter 
(longer) baseline there should be lesser (larger) NSI effect. If the experimental data is available over certain 
range of neutrino energy from the shorter and longer baseline experiments then only the multi-parameter ht 
will help in disentangling effects due to different unknown parameters and hence improving bounds on NSI 
parameters or discovering it. 


In the following we present three subsections based on three different sources of neutrino flux : Superbeam 
in HI A, Neutrino Eactory in HIB and monoenergetic neutrino beam in HI C. Particularly for Superbeam and 
Neutrino factory, we have studied the NSI effect in CP violation discovery in the leptonic sector in neutrino 
oscillation experiments. Monoenergetic neutrino beam requires shorter baseline due to technical reasons as 
discussed later. Due to that the NSI effect in the monoenergetic neutrino beam studies is insignificant and 
we will be discussing only CP violation discovery reach in absence of NSI. Two different nuclei have been 
considered for monoenergetic beam and to study CP violation. 
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A. Superbeam 


We discuss below the experimental set-up and systematic errors for superbeam. We have presented our 
results for discovery reach of CP violation due to 6 and other off-diagonal NSl phases (f>ij for different 
experimental set-up with different baselines. We have shown also the effect of the absolute values of different 
£a /3 in CP fraction. 


1. Experimental set-ups and systematic errors 


We consider three experimental set-ups: (a) A Superbeam set-up originating in CERN and reaching a 
500 kt Water Cherenkov detector [79] placed at a distance of 130 km at Frejus (France), (b) A Superbeam 
set-up originating in Tokai and reaching a 500 kt Water Cherenkov detector [13, 80, 81] placed at a distance 
of 295 km at Kamiokande (Japan) and (c) A Superbeam set-up which originates in CERN and reaches a 100 
kt Liquid Argon detector placed at a distance of 2300 km at Pyhasalmi (Finland). The observable channels 
that we have considered are and —)■ 

The flux considered for set-up (a) has mean energy ^ 0.3 GeV, which are computed for 3.5 GeV protons 
and 10^^ protons on target per year. For our analysis the beam power has been considered of about 4 MW per 
year and the time period has been taken to be 2 years for neutrinos and 4 years for anti-neutrinos. We consider 
the same flux as in [14, 82]. In the case of set-up (a) the efficiencies for the signal and background are included 
in the migration matrices based on [79] except for the channels disappearance, disappearance and 
(NC) which are 64%, 81% and 11.7% efficiencies respectively. We have considered systematic uncertainties 
of 2% on signal and background channels. 

The experimental set-up (b) has been considered as given in GLoBES [83] and the detector specification 
is mentioned in table III: 


Target Power 

4 Mega-Watt 

Fiducial Mass 

500 kt 

Data Taking 

4 yrs v and 4 yrs V 

Baseline 

295 km 

Energy Resolution 

0.085 


TABLE III: Detector characteristics for Hyper-Kamiokande. 

The flux considered for set-up (c) has mean energy ^ 5 GeV, which are computed for 50 GeV protons 
and 3 x 10^^ protons on target per year. For our analysis the beam power is about 0.8 MW per year and 
the time period has been taken to be 5 years each for neutrinos and anti-neutrinos. We consider the same 
flux as in [18]. The detector characteristics for set-up (c) is given in table IV [84]. The correlation between 
the visible energy of background NG events and the neutrino energy is implemented by migration matrices 
which has been provided by L. Whitehead [85]. 

In table V for <5 = 0 we have shown the expected number of events for three experimental set-ups, for no 
NSIs and also for different real NSIs , each of which equal to a (one NSI at a time). We have considered 
the central values of various parameters as shown in table II and have used PREM proHle [75] for matter 
densities. One can see that there is significant effect of and Eer on the number of —)■ events. 
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Signal Studies 

Ve CC appearance Studies 

Vfj^ CC Disappearance Studies 

Signal efficiency 

80% 

85% 

NC mis-identification rate (Background) 

1% 

0.5% 

CC mis-identification rate (Background) 

1% 

0% 

Signal Normalization error 

5% 

10% 

Background Normalization Error 

15% 

20% 

Neutrino Energy Resolution 

Ve CC energy resolution 

0.15y/A(GeV) 

CC energy resolution 

0.2i/E(GeV) 

scale uncertainty 

2% 

scale uncertainty 

0.01% 


TABLE IV: Detector characteristics for set-up (c). 



CERN-Frejus 
(130 km) 

T2HK 
(295 km) 

CERN-Pyhasalmi 
(2300 km) 

SM 

4533 

6098 

421 

^ee 

4536 

6110 

426 


4634 

6272 

502 

^er 

4543 

6143 

470 

^ flT 

4529 

6081 

414 


4533 

6099 

425 

Srr 

4530 

6085 

413 


TABLE V: The number of —)■ events for CERN-Erejus (130 km), T2HK (295 km) and 

CERN-Pyhasalmi (2300 km) for SM and for different NSI (one at a time each of which equals to a). 


Possibilities of significant effect of these NSIs can be seen in the expression of oscillation probabilities in 
Eq.(lO). 


2. Results 


1000 



CERN - 

FREJUS 

-- 

T2HK 


■■■ 

CERN - 

Pyhasalmi 


100 



6 (true) 


FIG. 1: Discovery reach of CP violation due to 5 for three different experimental set-up : CERN-Frejus, T2HK 
and CERN-Pyhasalmi set-up considering only SM interactions for normal hierarchy. 
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FIG. 2: F(5) for three different experimental set-up: CERN-Frejus, T2HK and CERN-Pyhasalmi set-up at So- 

considering real NSIs See, and Ett- 


Here "we first discuss the CP violation discovery reach "without any NSI for the three experimental set-ups 
(a), (b) and (c). Then we make a comparative study between no NSI and with real NSI. We also consider 
off-diagonal NSIs with phases. The Figures are given only for normal hierarchy (true). However, there is no 
significant departures in the results for inverted hierarchy from that with normal hierarchy. 

In Figure 1 we have shown the values versus i5(true) for SM interactions of neutrinos with matter 
from which the discovery reach can be obtained at different confidence levels. We have fixed (5(test) to 0 
and TT and have marginalized over hierarchy for every 5(true) value. To calculate CP fraction from Figure 
1 at 3a one is required to find the fraction of the total allowed range of 6 for which values are above 
the solid horizontal line in Figure 1. Particularly at 3cr confidence level for CERN-Frejus set-up the CP 
violation could be discovered with F{6) of around 0.58 of the possible <5 values for normal hierarchy (true 
value) whereas for T2HK set-up these values are around 0.67 and for CERN-Pyhasalmi set-up these values 
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0.8 - 

« 0.6 - 

'©■ 

fa 0.4 - 

0.2 - 

0.0 — 
0.0001 


ej,i(true) 

FIG. 3: F{(j>sii), F{(j)sT) and F{(f>i_n-) (one NSl at a time) for 5cp = 0 at Scr confidence level for three different 

experimental set-up. 

are about 0.21. For longer baseline for normal hierarchy the discovery reach is little bit better than the 
inverted hierarchy (not shown in the Figure). The discovery reach for longer baseline of 2300 km was shown 
earlier by Coloma et al [18]. So with only SM interactions of neutrinos with matter the short baseline like 
T2HK set-up seems to be better for good discovery reach of CP violation. This was observed earlier by 
different authors [18, 19, 24-26]. 

a. Possibilities of Discovery of CP violation due to 6 for real NSI: If we compare the F{5) 
obtained from Figure 1 at 3a with that shown in Figure 2 at 3 ct we find that for smaller values of NSIs 
satisfying model dependent bound for short baseline of CERN-Frejus and T2HK set-up , there is insignificant 
effect of NSI on the the discovery reach of CP violation due to d. However, for CERN-Pyhasalmi set-up 
even for smaller NSI values one finds the effect on CP violation discovery reach. One interesting feature is 
found from these figures which is that if one does not observe CP violation due to 6 in shorter baselines (say 
CERN-Erejus set-up) there is still a possibility of observing CP violation in longer baselines (see CERN- 
Pyhasalmi set-up with off-diagonal NSI like Se/x and Ser) which could be the signal of NSI with significant 
strength. In Eigure 2 one can see that except for other NSIs the F{S) changes considerably, with increase 
in the values of different NSIs for CERN-Pyhasalmi set-up. 

b. Possibilities of Discovery of CP violation due to complex NSI phases : There could be CP 
violation only due to off-diagonal NSI phases for which 5 = 0. Such cases have been shown in Eigures 3. 
We observe that for certain NSI value the CP fraction for longer baselines is more in comparison to that 
for the shorter baselines. With the increase of \eap\ there is increase in discovery of CP fraction in general. 
If CP violation is not observed in short baselines (say CERN-Frejus) but it is observed at relatively longer 
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ItepI = 3.8 x10 Iee,| = 0.25 Itp^l = 0.25 
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0 e/.i (true ) (/»et (true ) (true ) 


|£,^| = 3.8 XiO kerl = 0.25 |£pr| = 0.25 



4 ’.,, (true ) 0.1 (true ) 0 ., (true ) 

FIG. 4: Allowed region (unshaded) for total CP violation discovery reach F{S) for two different baselines T2HK 
(top panel) and CERN-Pyhasalmi set-up (bottom panel) considering NSI phases (Aer and 

baselines (say CERN-Pyhasalmi set-up) then it would signal the presence of NSI as can be seen from Figure 

3. 

c. Possibilities of Discovery of CP violation due to both 5 and complex NSI phases : In 

Figure 4 in the un-shaded region discovery of total CP violation could be possible at 3 (t confidence level. 
Here, we are assuming that nature has two sources of CP violation - one due to 5 phase and other due to 
one NSI phase. For this we ask whether there is (J — (pij) parameter space where both can be distinguishable 
from their CP conserving points as discussed earlier. From the Figure it is seen that for T2HK set-up there 
is more allowed region and as such better discovery reach in the 5 — (j)ij parameter space particularly for 
NSI - Sell in comparison to CERN-Pyhasalmi set-up . Particularly for NSI - Ser and the discovery 
reach of total CP violation is better for both longer and shorter baselines (like CERN-Pyhasalmi set-up 
and T2HK set-up). Simultaneous consideration of experimental data from both these set-ups could improve 
the discovery reach further because of presence of non-overlapping un-shaded discovery regions in the upper 
panel two right Figures and lower panel two right Figures corresponding to two different NSIs. 

From Figure 1 it is observed that for certain range of values of S (about 0°to 25°, 160° to 208°, 340°to 360°), 
the CP violation due to 6 could remain unobservable. However, even if nature has such values of S, from 
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Figure 4 one can see that total CP violation could be observed in such cases in presence of NSI for a wide 
range of values of different NSI phases. However, disentangling the sources of two types of CP violating 
phases (one from NSI phase and other from S) could be difficult. Here the total CP fraction can be thought 
of as the ratio of the un-shaded region divided by the total region covered by <5 — </)ij parameter space. As for 
example, in the upper panel extreme left Figure such total CP fraction is about 75%. For both the set-ups 
for SgT- such CP fraction is even more. For CERN-Pyhasalmi set-up for Sg^ the total CP fraction is not 
good. 


B. Neutrino factory 


We discuss below the experimental set-up and systematic errors for neutrino factory. We have presented 
our results for discovery reach of CP violation for different baselines of length 730 km (FNAL-Soudan),1290 
km (FNAL-Homestake) and 1500 km (FNAL-Henderson) due to S and other off-diagonal NSI phases 
We have shown the effect of the absolute values of different Sajs in that. For this analysis we have considered 
one NSI at a time. 


i. Experimental set-ups and systematic errors 


We have used a large magnetised iron neutrino detect or (MIND) [44] with a toroidal magnetic field having 
a mass of 100 KTon. MIND can also be described as an iron-scintillator calorimeter. This detector has 
the capability of excellent reconstruction and charge detection efficiency. In this section we have considered 
muons in a storage ring consisting of both and which decay with energies of 10 GeV. We consider 5 x 
10^^ stored muons. The golden channel (i^g —)■ oscillation channel) where the charged current interactions 
of the Vfj, produce muons of the opposite charge to those stored in the storage ring (generally known as wrong- 
sign muons), is the most promising channel to explore CP violation at a neutrino factory. The detector that 
we are considering in this work - MIND is optimized to exploit the golden channel oscillation as this detector 
has the capacity to easily identify signal i.e. a muon with a sign opposite to that in the muon storage ring. 
Different oscillation channels which have been considered as signals and backgrounds [44] in the analysis are 
shown in table VI. We have taken the migration matrices for the true and reconstructed neutrino energies 
as given in reference [44]. The signal and background efficiencies are taken into account in those matrices. 
We have taken into account the reconstruction of t contamination coming from the silver channel Vg — )■ v-g 
as background. We have considered systematic errors to be 1%. In this work we have considered a running 
time of 5 years for both and [i~ . 



Channel Name 



Signal 

Golden Channel 

Vg V^g 

Vg Vfj, 


Vg disappearance channel 

v^g Vg 

Vfg Vg 


Silver Channel 

Vg Vr 

Vg Vr 

Background 

v^g disappearance channel 

Vfl -t Vfl 

Vfg -t v^ 


Platinum Channel 

Vfl Vg 

Vfg Vg 


Dominant Channel 

Vfg -t Vr 

Vfg Vr 


TABLE VI: Different oscillation channels considered as signals and backgrounds in the analysis. 
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730 km 

1290 km 

1500 km 

SM 

248491 

222097 

210819 

^ee 

248715 

222644 

211527 


316049 

283709 

269412 

^er 

251331 

229148 

220134 


248049 

220985 

209394 


248682 

222667 

211597 

£rr 

248072 

220964 

209313 


TABLE VII: Number of Vf, —>■ Vy, events for 730 km (FNAL-Soudan),1290 km (FNAL-Homestake) and 1500 
km (FNAL-Henderson) for SM and for different NSI (one at a time each of which equals to a). 


In table VII for 5 = 0 we have shown the expected number of events for three baselines in neutrino factory, 
for no NSI and also for different real NSI , each of which equal to a (one NSI at a time). We have considered 
the central values of various parameters as shown in table II and have used PREM profile [75] for matter 
densities. Like superbeam case here also for neutrino factory one can see that there is significant effect of 
£f^y and £eT on the number of v^. —>■ Vy events. Possibilities of significant effect of these NSIs can be seen in 
the expression of oscillation probabilities z/g —> (which can be obtained using Eq.(13) in Eq.(lO)). 


2. Results 
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FIG. 5: F{S) with only SM interactions of neutrinos with matter for baselines of different lengths and different 

muon energies Ey at 5 ct confidence level. 


The optimization of CP violation discovery reach has been done earlier [29, 30] for different baselines and 
different parent muon energy when only SM interactions of neutrinos with matter during propagation is 
present. We have re-analysed this optimization based on the updated MIND detector characteristics and the 
migration matrices as given in [44]. We present the result for neutrino factory at 5cr confidence level unlike 
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superbeam case where we have given the results at 3a confidence level. This is because MIND detector for 
neutrino factory is found to give quite good CP violation discovery reach. It is found that at 5cr confidence 
level the CP fraction (shown as F{S) in Figure 5) of about (0.8 < F{S) < 0.9 ) is possible for baselines 
ranging from about 300 to 800 km for energies lesser than 5 GeV and for baselines ranging from about 700 
to 2000 km for energies 6-10 GeV respectively. Based on high CP fraction discovery potential as found in 
this Figure we have chosen 10 GeV muon energy and a few baselines which are : 730 km(FNAL-Soudan), 
1290 km (FNAL-Homestake) and 1500 km (FNAL-Henderson). 
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FIG. 6: F{S) versus real NSI (ep) at 5a confidence levels. 


Next we have tried to answer the question that had there been non-Standard interactions what could 
have been their effect on the CP violation discovery reach for such experimental set-ups. While taking into 
account NSI effect, for off-diagonal NSI we have also considered the effect of NSI phases over 5cp violation. 
In contrast to the results on superbeam we have presented results for neutrino factory at 5cr confidence level 
in Figures 6 and 7 because in most of the cases the discovery reach of CP violation is far better in neutrino 
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factory set-ups than those in superbeam set-ups particularly for NSIs ^ 1. However, for comparison of the 
results of neutrino factory set-ups with that of superbeam set-ups we have considered the 3 ct results for both 
type of set-ups. For brevity we have avoided presenting the results at Scr corresponding to Figures 6 and 7. 

a. Possibilities of Discovery of CP violation due to 5 for real NSI: In Figure 6, we have studied 
F{5) in the presence of real NSI (NSI phases have been chosen to be zero) for different baselines of length 730 
km, 1290 km and 1500 km at 5 (t confidence level. Here, in plotting the Figures we have considered the NSI 
values upto the model independent bounds as shown in table I. For lower values of NSI there is essentially 
negligible effect on discovery reach of CP violation which is seen in the Figure as horizontal straight line. 
This corresponds to F{5) due to SM interactions only which can be verified from Figure 5 at muon energy of 
10 GeV for the appropriate baseline. For NSI with model dependent bound there is insignificant effect 
on F{5) whereas for other NSIs there are some effect on F{5). Except ee^, Eer and for other NSIs there 
is considerable decrease in discovery reach of CP violation with the increase of NSI values. Even if there is 
CP violating phase 6 but as one can see that in presence of higher allowed values of NSI - e^-t there could 
be no discovery of Jcp violation. 

The results are shown in Figure 6 at 5a confidence level. However, at 3cr also the decreasing features of 
Eee and Ett ^ 1 is present and in comparison to Figure 2 it is found that for large NSI values of Egg, £tt ^ 1 
superbeam T2HK set-up is found to be better for the discovery reach of CP violation as compared to neutrino 
factory set-ups. But for lower values of NSIs particularly £ < 10“^ the neutrino factory set-ups are found to 
be much better for the discovery of CP violation at ia and even at 5a the CP fraction is more in neutrino 
factory set-ups. 
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FIG. 7: F{(j)gr) and F{(j>i_i,T) for 5cp = 0 at 5 ct confidence level for three baselines. 
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FIG. 8: 5 versus phase {(juj) considering the value of NSIs(eij) at the uppermost value with model dependent 
bounds. Unshaded regions correspond to CP violation discovery reach. 


b. Possibilities of Discovery of CP violation due to complex NSI phases : In Figure 7 we have 
addressed the question of what could be the CP fraction for discovery of CP violation if Dirac phase <5 is 
absent in PMNS mixing matrix and CP violation comes from purely NSI phases. We observe that for certain 
NSI value the CP fraction for longer baselines could be more in comparison to that for the shorter baselines. 
With the increase of |£a/ 3 | there is increase in discovery of CP fraction in general. If CP violation is not 
observed in short baselines (say 730 km) but it is observed at relatively longer baselines (say 1500 km) then it 
would signal the presence of NSI as can be seen from our Figures although for it is difficult to distinguish 
the effect between short and long baselines considered by us. Comparing 3 ct result of neutrino factory (not 
shown here) with superbeam result of Figure 3 it is found that in the presence of £e^ the Neutrino factory 
set-up gives better CP violation discovery reach whereas for e^r the superbeam CERN-Pyhasalmi set-up is 
better . For Ser the neutrino factory set-up is slightly better than superbeam set-up. 
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c. Possibilities of Discovery of CP violation due to both 5 and complex NSI phases : In 

Figures 8 we have considered the case for 730 km baseline where the CP violation might come from 5cp as 
well as from NSI phase ipap- For 1290 km and 1500 km baselines the essential feature is almost same (not 
shown in the Figure). Here we have chosen the uppermost value of NSI with model dependent bound. The 
total CP violation due to S and pij could be observable even for those S which could remain unobservable 
in absence of NSI phases for some ranges of S values. In fact, the CP violation discovery reach could be 
possible for the whole range of allowed S values for some specific values of pij as can be seen from Figures. 
But it would be difficult to conclude whether the observed CP violation is due to 6 or any NSI phase. As 
compared to Figure 4 due to different superbeam set-ups, we find better discovery reach of CP violation for 
neutrino factory set-up for 730 km baseline in Figures 8 . As for example, in the extreme left Figure for Sefi 
the total CP fraction is around 94% for neutrino factory in contrast to 75% for the T2HK superbeam set-up 
as shown in Figures 4. 

Comparing Figures 8 and 4 we find that for if we combine experimental data from superbeam set-ups 
and neutrino factory set-ups, there are scopes to improve the total CP violation discovery reach further. 


C. monoenergetic neutrino beam 


We discuss below the possible discovery reach of CP violation for two different nuclei, 5 o°Sn and ^®^Yb 
which is considered for the electron capture experiment. We also discuss our procedure of energy resolution 
at the detector for numerical simulation. Earlier in section II we have shown how —)■ oscillation 
probability depends on 5 for shorter baseline (A ~ a). Based on the oscillation probability and its’ variation 
with respect to 6 we discuss the procedure for choosing suitable boost factor 7 numerically for specific 
baseline and specific nuclei considered for the neutrino beam. We mention below four experimental set-ups 
that we have considered for analysis of discovery reach oi CP violation with monoenergetic neutrino beam. 


1. Suitable boost factor, neutrino energy from Ve —> oscillation probability 


The most suitable candidate for producing neutrino beam from electron capture process would be the one 
with a low Q value and high boost factor, 7 [48]. Also it would be preferable if the nuclei has a short half life. 
The reason for these is as follows. We need neutrino energy around the peak of the oscillation probability 
where variation due to 5 is significant and as such « (2n -b l)f. Considering E = 2jQ, it follows 

that 7 = as for example, for the first oscillation peak. For sufficiently high 7 almost all neutrinos 

are expected to go through the detector. Then to satisfy the above condition we need to lower Q value. 
Then another condition is, ■jt < T where T is the time considered to perform the experiment so that all the 
nuclei decay and r is the half life of the nuclei. If 7 is increased then the half life r is required to be small. 
So the preferable factors considered in choosing the candidate for producing neutrino beam from electron 
capture process are - low Q value, small half life t and high 7 [48]. Although higher 7 needs technological 
advancement for the accelerator. 

The isotope, sJ^Sn, has Q = 267 keV in the rest frame and a half life of 4.11 h. As it has a low Q value 
so one may consider high 7 value. However, it has a longer half life as compared to other nuclei like 
^®^Yb, ^^®Yb, ^^®Dy [47, 50] whose half lifes are 18.5 seconds, 3.04 seconds, 261 seconds, 7.2 min and 
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3.1 min respectively. However, these nuclei have larger Q values of the order of 10^ keV. On the other hand, 
considering effective running time per year as 10^ second all the nuclei for isotope 5 o*^Sn will not decay. But 
for 7 = 500 or 320 (as considered in our analysis to obtain the suitable neutrino energy E resulting in high 
oscillation probability) the useful decays are respectively about 0.608 or 0.768 times the total number of 
IJ^Sn nuclei considered. So there is not much suppression in numbers of nuclei. Hence although IJ^Sn has 
a larger half life, due to its lower Q value there is scope to consider higher 7 for CERN-Frejus or 250 km 
baselines. For these reasons we have preferred isotope, sJ^Sn in comparison to other nuclei. However, there 
is recent study on Hnding suitable candidate nuclei for electron capture process and it has been found that 
^^^Yb, ^^®Yb nuclei have dominant electron capture decay to one level. Particularly, ^®^Yb has been 
found to be most suitable one [50]. For this reason, apart from nuclei 5 o°Sn we shall consider ^®^Yb also for 
our analysis. However, as Q value (5435 keV) for ^®^Yb is higher, corresponding 7 value for such nuclei are 
supposed to be small. 

The neutrino beam produced from electron capture process is boosted with certain boost factor, 7 . The 
boosted neutrino beam produced from such process hits the detector at a baseline of length L at a radial 
distance R from the beam axis and the energy, E of this beam in rest frame of the detector , i.e., in laboratory 
frame is given by: 


E{R) 


QL_ P 

7 [ Vl + (R/L)2. 


1 + (7R/L)2 


( 22 ) 


where R is the radial distance at the detector from the beam axis. At beam center, R = 0. From the above 
equation ( 22 ) the energy window considered for the analysis which is constrained by the size of the detector 
is given by: 


27Q 

1 + {iR-max/LY 


<E< 27 Q 


(23) 


From equation (23) we can see that once the baseline length L and 7 is Hxed the energy window gets fixed. 
However, even considering radius of the detector Rmax = 100 m the energy window is very small as can be 
seen from Figure 9. 


One can see from equation Eq.(22) that it is possible to tell precisely the energy from the R value of 
the Cherenkov ring at the Cherenkov detector instead of measuring directly the neutrino energy. So there 
is scope to get good energy resolution by measuring position if the vertex resolution is good. The <j{E) 
function corresponding to energy resolution function (as used in running GLoBES [76]) in terms of vertex 
measurement uncertainty cr(i?) can be written as: 


(j{E) 


_ QRp 

+ (1- 



— cr{R) 

7 


(24) 
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where j3 is defined as 




- 1 

7 


(25) 


Vertex measurement uncertainty for electron (muon) identification at Super-K is around 30 (20) cm [86]. To 
estimate the energy resolution using position measurement one may consider a{R) ^ 30 cm provided that the 
beam spreading is negligible ( lesser than about 1/rrad) [45] which is difficult to achieve experimentally. 
If we take into account the beam divergence about 10 fxrad (which is almost one order larger than that 
considered in references [87, 88]), one may consider larger a{R) ^ 130 cm particularly for baseline of 130 
km. For baseline like 250 km it would be more but we have considered same u{R) which means the beam 
divergence has been assumed to be lesser than about 5 fxrad for the analysis for baseline with length 250 
km. 


In this work GLoBES[76] has been used for doing the simulations. In order to use this software, the radial 
binning is replaced by binning in energy and the bins are not equidistant. If we divide R^ax into k bins the 
edges of the bins are given as: 


= R^ax - (t - l)Ai?2 


(26) 


with 


Ai?2 


d 2 

-^mi 

~k 


(27) 


We consider Rf > 


so that in GLoBES the respective energy bins are in the correct order as given below 


E{Rf) < E{Rl,) 


(28) 


where E is the neutrino energy in the lab frame. 


The number of events per bin i and channel c (different channels mentioned later in this section) is given 
by: 


pT pEi+AEil2 poo 

dE' dE(j){E)P^{L,E)a^{E)e^{E')R^{E,E') (29) 

T JEi-\Ei/2 Jo 

where Nnorm is the normalization factor for using GLoBES and is related to length of the baseline, area 
and energy binning related to flux , number of target nuclei per unit target mass and number of nuclei 
decaying, is the signal efficiency in the respective channel, P‘^{L,E) is the neutrino oscillation probability 
in particular channel, cr‘^(E) is total cross section for particular flavor of neutrinos and particular interaction 
corresponding to particular channel. R{E, E') is the energy resolution function of the detector where E' 
is the reconstructed neutrino energy. (j){E) has been calculated from the angular neutrino flux ^{E) as 



defined below : 
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_ ^decays f ^on\ 

^ “ Att \q) ^ ^ 

where Ndecays is number of nuclei actually decaying per year. The detailed derivation of these expressions 
can be found in [45]. Considering equation (23) and ( 30) one can see that with increase in 7 value the 
angular flux increases. 



Energy (GeV) 



Energy (GeV) 


FIG. 9: Probability P(i^e —>■ Vy) vs neutrino energy E for two different nuclei so^Sn and ^®^Yb and the corresponding 
energy window satisfying equation (23) for different 7 values. 


We have plotted numerically the probability P{ve t'/r) with respect to energy for two different baselines 
of length 130 km (CERN-Frejus) and 250 km for two different nuclei 5 o°Sn and ^®^Yb. We have considered 
normal hierarchy in plotting Figure 9. For all plots in Figure 9 , <5 has been varied over its’ entire allowed 
range (0 to 2 tt) resulting in the shaded region in each plots showing the significant variation of probability 
at particular neutrino energies. Corresponding to each of the nuclei (whose Q values are fixed) we have 
considered appropriate 7 value so that the corresponding energy window (as mentioned in (23)) overlaps 
with the shaded region near the oscillation peaks having significant variation of probability due to 5 variation. 
In choosing 7 , one also has to keep in mind that the neutrino energy is not too low as otherwise flux will be 
much lesser. The energy window has been shown by the shaded vertical strips. For our suitable choice of 7 
value, the energy window is larger for 130 km baseline and relatively smaller for 250 km baseline for both 
the nuclei. Also the energy window for Ip^Sn is larger than ^®^Yb. 

For finding 5 we shall prefer the maximum variation of the probability with 5 which will occur for neutrino 
energy satisfying —« (2n + l) 7 r /2 where n is an integer. This has been shown in Figure 9 in which 
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the oscillation probability has been evaluated numerically considering the evolution of neutrino flavor states. 
However, the energy also depends on the Q value of the corresponding nuclei. So we have considered the case 
of two nuclei separately. As for example, for 5 o°Sn nuclei, for baseline of length 130 km we have considered 
first oscillation maximum and for baseline of length 250 km we have considered second oscillation maximum 
where the variation of the probability due to 6 is significant. For ^^^Yb, for both the baselines we have 
considered first oscillation peak. In considering the suitable peak in the oscillation probability one has to 
keep in mind that the neutrino flux increases with as shown in Eq.(30) and so after doing the numerical 
analysis one can decide which energy out of various energies near various peaks are suitable. However, one 
also have to think about feasible boost factor 7 . For that reason although the first oscillation peak for ^J^Sn 
nuclei corresponds to higher energy with respect to 2 nd oscillation peak but as it also requires higher boost 
factor 7 around 900, we have considered the neutrino energy near the second oscillation peak for 5 o*^Sn nuclei 
for baseline of length 250 km. Depending on the energy chosen near a peak one can appropriately choose the 
boost factor 7 on which the neutrino energy window as shown in equation (23) as well as Vg flux as shown 
in equation (30) depend. This has been illustrated in Figure 9. 


2. Experimental set-ups and systematic errors 


For doing the analysis we choose four different set-ups: 
set-up(a): The length of the baseline is taken to be 130 km (CERN-Frejus baseline) and the boost factor 
7 to be 500 for nuclei sJ^Sn. 

set-up(b): The length of the baseline is taken to be 250 km and the boost factor 7 to be 320 for nuclei 
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set-up(c): The length of the baseline is taken to be 130 km (CERN-Frejus baseline) and the boost factor 
7 to be 26 for nuclei ^^^Yb. 

set-up(d): The length of the baseline is taken to be 250 km and the boost factor 7 to be 40 for nuclei ^^^Yb. 

We consider a Water Cherenkov detector of fiducial mass 500 kt. Following [89], the signal efficiency is 
considered to be 0.55 for appearance channel. Background rejection factor coming from neutral current 
events is considered to be 10“^ for appearance channel. Signal error of 2.5% and background error of 5% 
has been considered. For quasi-elastic appearance and Ve disappearance we have followed signal efficiency 
and error as given in reference [89]. We have considered the neutrino energy resolution as discussed earlier in 
(24) which can be obtained from vertex resolution after taking into account beam spreading. The neutrino 
energy is known from Eq.(22) and the energy width considered by us is obtained from Eq.(23) by considering 
the radius of the detector Rmax = 100m. We assume 10^® electron capture decays per year and the running 
time of 10 years for accumulating data. 

However, depending on the half life of Sn (4.11 hrs), the number of useful decays per effective year (10^ 
seconds ) considered are about 0.608 x 10^® with boost factor (7 = 500) for 130 km baseline and 0.768 x 10^® 
with boost factor (7 = 320) for 250 km baseline. Also, depending on the half life of ^®^Yb (3.04 seconds) 
and the boost factor 7 = 26 or 7 = 40, for baselines 130 km or 250 km respectively, the number of useful 
decays per effective year ( 10 ^ seconds ) considered are almost equal to the total number of nuclei i.e, 10 ^® 
as half life is much smaller than sJ^Sn. It is possible to achieve 7 about 480 at upgraded SPS facility at 
CERN [90-92] and 7 > 1000 for LHC based design [93]. We have considered six energy bins keeping in mind 
the available energy window for different set-ups and the corresponding energy resolution in equation (24). 
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In considering the energy resolution we have taken into account beam spreading. For that the the energy 
resolution considered by us is bad in comparison to the energy resolution considered in reference [45] and we 
have to consider much lesser number of energy bins. 

As the number of events corresponding to all set-ups are quite large, any background due to atmospheric 
neutrinos are expected to be quite small and we have not considered such background in our analysis. 


3. Possibilities of Discovery of CP violation 


Here we discuss the discovery of CP violation for four different experimental set-ups (a-d) mentioned earlier 
for monoenergetic neutrino beam. In presenting our analysis we have followed the numerical procedure as 
discussed at the beginning of section III. We have considered the true hierarchy as normal hierarchy. 
However, we have considered the uncertainty in the hierarchy of neutrino masses in the test values as it 
is not known at present. We have also considered the uncertainties in the other oscillation parameters as 
mentioned in table H. For finding CP violation we have fixed (I(test) at CP conserving 5 values (0,7r). 

In Figure 10, versus 6 (true) has been plotted to show the discovery reach of the CP violation for two 
different set-ups - set-up(a) and set-up(b) for |o°Sn nuclei. We find that the discovery of CP violation for 
set-up(a) & (b) could be found with F{5) for about 51% and 49% respectively of the possible 6 values at 3 ct 
confidence level. In Figure 11, versus 5 (true) has been plotted to show the discovery reach of the CP 
violation for two different set-ups - set-up(c) and set-up(d) for ^^^Yb nuclei. We find that the discovery of 
CP violation for set-up(c) and set-up (d) are not good and could not be found even at Icr confidence level. 
For longer than 250 km baselines we have not presented any plots for CP violation discovery reach. It seems 
one of the basic problem for longer baselines will be relatively bad energy resolution because we are trying 
to use vertex resolution for getting energy resolution but there is beam spreading and as such over longer 
baseline beam spreading will make the energy resolution poorer. 



FIG. 10: versus 5 (true) for two experimental set-ups (a) & (b) with nuclei 5 o°Sn for true normal hierarchy. 
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5 (true) 5 (true) 


FIG. 11: Ax^ versus (5(true) for two experimental set-ups (c) & (d) with nuclei 


i52Yb foj. true normal hierarchy. 


In our analysis we have chosen neutrino energy near the oscillation peak (as shown in Figure 9) which is 
more S sensitive region. This consideration improves the CP violation discovery reach. However, in [45] (as 
shown in Figure 7 of that paper) CP violation discovery reach has been shown to be about 81% of the possible 
S values for their set-up II for 250 km baseline for the presently known 0 i 3 value. The discovery reach seems 
to be not so good in our case although we have chosen appropriately the neutrino energy where probaility 
of oscillation is <5 sensitive. It is because in our analysis we have considered more realistic value of 7 (which 
is lower than that considered in [45] but could be achievable at present (keeping in mind the possible SPS 
upgrade at CERN). We have also taken into account the beam spreading in our present analysis in estimating 
effective energy resolution and have taken only a few neutrino energy bins for the analysis unlike [45] as it 
is important to consider the size of the energy bins larger in comparison to the level of energy resolution. 
Our results shows that the nuclei 5 o°Sn is more suitable than the nuclei ^^^Yb so far CP violation discovery 
reach is concerned. 


IV. CONCLUSION 


We have studied the possible CP violation discovery reach due to Dirac phase S in the leptonic sector 
through neutrino oscillation experiments with superbeam, neutrino factory and mono-energetic neutrino 
beam from electron capture process. Particularly for superbeam and neutrino factory we have studied 
the NSI effect on the CP violation discovery reach. However, for mono-energetic neutrino beam we have 
considered shorter baseline as required technically for lower boost factor 7 . For that NSI effect seems to be 
insignificant and hence we have studied discovery of CP violation due to S without different NSI in this case. 

For short baseline with CERN-Frejus and T2HK set-up in case of superbeam for different NSI satisfying 
model dependent bound we find that there is insignificant effect due to real NSI on the the discovery reach 
of CP violation due to S. For longer baseline with CERN-Pyhasalmi set-up such effects are significant. In 
case of neutrino factory, the baselines considered are slightly longer than those for superbeam. In this case, 
except other NSIs have significant effect on F(S). 

As P(S) is never found to be 1, obviously for some values of <5, CP violation may not be observed at 
certain confidence level for the kind of detectors we have considered (see Figure 1 ). Even if one observes CP 
violation due to S in shorter baselines there is a possibility of observing no CP violation in presence of real 
NSI. This feature can be seen in Figure 2 for superbeam in presence of real NSIs - e^e, £tt and in Figure 6 
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in presence of real NSI Err for neutrino factory. For smaller values of NSI, the CP violation discovery reach 
is much better for neutrino factory set-ups than those for superbeam set-ups. However, if NSIs like See, Ser 
and Ett are of significant strength ^0.1 then CP violation discovery reach at neutrino factory could be very 
bad in comparison to particularly T2HK set-up in superbeam. One interesting feature is found from Figure 
2. If one does not observe CP violation due to 5 in shorter baselines (say CERN-Frejus set-up) there is still a 
possibility of observing CP violation in longer baselines (see CERN-Pyhasalmi set-up with off-diagonal NSI 
like Eefi and e^t) which could be the signal of NSI with significant strength. Using short and long baseline 
one could conclusively tell about CP violation due to 5 and about NSI under such situation. 

NSI - Eefi, Eer and E could be complex. We have considered the corresponding phases (pefi, 4>eT and (p^r 
respectively in the analysis of the discovery reach of total CP violation for both in absence and in presence 
of 5. Even in absence of 6 one may observe CP violation due to the presence of NSI phases. The possibility 
of observing CP violation under such scenario is relatively better in general for longer baselines for both 
superbeam (see Figure 3 ) and neutrino factory (see Figure 7) provided that absolute values of NSIs are 
known. For Ee^ there is better CP violation discovery reach in neutrino factory set-up. For Est neutrino 
factory set-up with longer baseline is slightly better than CERN-Pyhaslmi set-up in superbeam at 3 (j. For 
Efir there is better CP violation discovery reach for CERN-Pyhaslmi set-up in superbeam. 

Assuming that there are two sources oi CP violation simultaneously existing in nature- one due to 5 and 
the other due to say one of the NSI phases (pij, the total CP violation discovery reach has been shown as 
unshaded region in the pij — 5 plane in Figure 4 for superbeam and Figure 8 for neutrino factory. From 
the Figures it is seen that for NSI - E^fi for T2HK set-up there is more allowed region and as such better 
discovery reach is possible in the 5 — pij parameter space particularly in comparison to CERN-Pyhasalmi 
set-up in superbeam. For NSI - Eer and the discovery reach of total CP violation is better for both 
longer and shorter baselines (like CERN-Pyhasalmi set-up and T2HK set-up). As compared to Figure 4 
due to different superbeam set-ups, we find significantly better discovery reach of CP violation for neutrino 
factory set-up for 730 km baseline in Figures 8. Comparing Figures 8 and 4 we find that for (per if we combine 
experimental data from superbeam set-ups and neutrino factory set-ups, there are scopes to improve the 
total CP violation discovery reach further. However, it would be difficult to disentangle the observed CP 
violation coming due to both 6 and NSI phase. In presence of some ranges of off-diagonal NSI phase values 
pefi and pfj^T (see upper panel of Figure 4 for superbeam and Figure 8 for neutrino factory) there is possibility 
of discovering total CP violation for any possible 5cp value at 3 (t. 

The CP violation due to 5 could remain unobservable with present and near future experimental facilities 
in superbeam and neutrino factory for certain range of values of 5 (as for example for superbeam with T2HK 
set-up in Figure I for about 0°to 25°, I60°to 208°, 340°to 360° at 3 cr confidence level). However, in presence 
of NSIs (with or without phases) the CP violation due to 5 or the total CP violation due to 6 and NSI phase 
could be observed even for such values of 5 as can be seen in various Figures shown in Section HI-A and B. 

Basic strategy to find CP violation in the leptonic sector in presence of NSI for superbeam set-up may be 
to consider both shorter baseline (say T2HK set-up) as well as one longer baseline (say CERN-Pyhasalmi 
set-up) because of their complementary nature with respect to the discovery reach oi CP violation. If 
NSI values are smaller (^ O.I) and real then the CP violation discovery in neutrino factory set-up with 
MIND detector seems significantly better than superbeam set-up. Eor complex NSIs, consideration of both 
superbeam set-up and neutrino factory set-up could give better CP violation discovery reach. 

We have discussed only the possibility of discovering CP violation. However, experimentally disentangling 
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CP violation coming from <5 present in PMNS matrix, the NSI phases and the absolute value of NSI could be 
very difficult. One may note that for longer baselines due to higher matter density NSI effect could be more. 
Also the effect due to d and that due to NSI on oscillation probability in particular channel varies differently 
with neutrino energy. Only if the experimental data is available over certain range of neutrino energy from 
the shorter and longer baseline experiments then only the multi-parameter ht could help in disentangling 
effects due to different unknown parameters. 

For monoenergetic neutrino beam sources we have considered two different nuclei - one Vg source is from 
electron capture decays of IJ^Sn isotopes and the other Vg source is from electron capture decays of isotopes 
^®^Yb. For each case we have considered two baselines 130 km and 250 km. Among experimental set-ups 
(a-d) in IIIC we find that the set-up (a) with sJ^Sn isotope and 130 km baseline is found to be the most 
suitable set-up for discovering CP violation with F{5) about 51 % at 3a confidence level. 

When one considers technical issues involved in the accelerator and running the ions through vacuum tube, 
isotopes ^®^Yb is better candidate than sJ^Sn isotopes because of much lesser half life. ^^^Yb is also better 
because of the dominant electron capture decay to one energy level. However, as can be seen from Figures 
10 and 11, the discovery reach of CP violation is found to be better for sJ^Sn isotopes. Out of different 
baselines for 5 Q°Sn nuclei, we find slightly better discovery reach for shorter baseline of 130 km with 7 = 500. 

Building up of such monoenergetic neutrino beam facilities will require some technological development 
and the implementation of it might take some time [94]. The existing CERN accelerator complex could be 
used to study such facility. However, the monoenergetic neutrino flux require a very large number of ions 
to be stored in the decay ring. It is difficult to control the beam at high intensities because of space charge 
detuning, intra beam scattering and vacuum loss. With SPS upgrade it could be possible to accelerate the 
ions to 7 = 480 but accelerating above that seems difficult [90-92]. Depending on the half life of sJ^Sn 
we have reduced the total number of useful decays of the ion per effective year from 10 ^® but the value 
considered is still extreme because of the requirement of acceleration and storage of the partially charged 
ion. For improving this the vacuum conditions in SPS would be required to be upgraded. It requires more 
study on such beam facility. With technological improvement if it is possible to consider monoenergetic 
beam with 7 > 1000 [93], then the CP violation discovery reach will improve further than what has been 
presented in this work. 
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